A Simple Method to Make Asymptotic Series of Feynman Diagrams Converge 
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We show that for two non-trivial Xcj)"^ problems (the anharmonic oscillator and the Landau- 
Ginzburg hierarchical model), improved perturbative series can be obtained by cutting oflt the large 
field contributions. The modified series converge to values exponentially close to the exact ones. 
For A larger than some critical value, the method outperforms Fade's approximants and Borel 
summations. The method can also be used for series which are not Borel summable such as the 
double-well potential series. We show that semi-classical methods can be used to calculate the 
modified Feynman rules, estimate the error and optimize the field cutoff. 
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Perturbative series associated with Feynman diagrams 
are commonly used in particle physics, solid state physics, 
optics and chemistry 1 1] . One remarkable success of this 
method is the prediction of the values of the anomalous 
magnetic moments of the electron and the muon with 
an incredible accuracy. Perturbative methods are also 
used to perform precision tests of the standard model of 
electro- weak and strong interactions Q. Despite these 
successes, it has been known for a long time [^,|| that 
series calculated from Feynman diagrams are not conver- 
gent but asymptotic. In other words, the range of validity 
shrinks with the order. One can improve this situation by 
using Pade approximants , either on the original series 
or a Borel sum Q of the series, if meaningful. However, 
even in the cases where the convergence of these alter- 
nate procedures can be proven, the convergence is very 
slow when the coupling is too large. In addition, for short 
series, it is difficult to estimate the error and to choose 
the best approximants. 

In this Letter, we construct improved perturbative se- 
ries which converge to values which are exponentially 
close to the exact ones. The error can be estimated 
analytically. The method can be applied on the lattice 
and in the continuum and works well when the meth- 
ods mentioned above are inefficient or not applicable. 
The method is a perturbative version of recent numerical 
calculations performed for various Xcf)^ models, namely 
the anharmonic oscillator |Q and the Landau-Ginzburg 
model in the hierarchical approximation In these 
calculations, we were led to introduce large (j) cutoffs and 
realized that as A increases, the field cutoff can be de- 
creased without affecting the accuracy of the result. The 
calculations presented here are perturbative series cal- 
culated with a large field cutoff. We only consider \4>'^ 
problems, however the procedure should extend to any 
kind of model where large field configurations are sup- 
pressed at positive coupling. This is in agreement with 
the general argument linking the large field configu- 
rations to the impossibility of applying Lebesgue domi- 
nated convergence to the path integral expression of the 
perturbative series. 

In order to give an idea about the efficiency of exist- 
ing methods, we consider the well-known example of the 



ground state energy of the anharmonic oscillator pO[ . 
The solid lines of Fig. ^ represent the number of signifi- 
cant digits obtained with perturbation theory for various 
A. As the order increases, the approximate lines rotate 
clockwise while moving left, forming an approximate en- 
velope. For a fixed coupling, there is an order of pertur- 
bation for which the error is minimized. On the other 
hand, the number of digits obtained with Pade approx- 
imants increases with the order. This is a consequence 
of Carleman's theorem which can be used to show pT| ] 
that diagonal sequences of Fade approximants converge 
to the ground state energy in an appropriately restricted 
domain of the complex plane. However, the convergence 
rate becomes slower as the coupling increases. This can 
be explained from the fact that, when the coupling 
increases, a [L/ L] approximant tends to a constant while 
the energy increases like A^/^. If Fade approximants are 
used for the Borel transform instead of the series, one 
obtains results qualitatively similar which are discussed 
later. 
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FIG. 1. Number of correct significant digits obtained with 
regular perturbation theory (solid lines) at order 1, 2, 3, 
15 and with Fade approximants [2/2], [3/3], .... [7/7] (dots) 
for the anharmonic oscillator, vs. logioA. The various orders 
can be identified from the explanations in the text. In all the 
graphs, the logarithms are in base 10. 

The method that we propose provides a systematic and 
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controllable improvement to regular perturbation and a 
better convergence than the Pade based methods in the 
right hand part of Fig. |l|. It can also be used in cases 
where the Borel sum has singularities on the positive real 
axis. In particle physics, the experimental error bars of 
precision measurements are shrinking and higher order 
terms of perturbative series are being calculated. Firmly 
established discrepancy or agreement between theory and 
experiment provides valuable information regarding the 
laws of nature at shorter distances. A common practice 
p2| to estimate at which order, for a given coupling, we 
reach the envelope illustrated in Fig. |], is to determine 
when the ratio of successive contributions reaches one. 
This method works quite accurately for the three exam- 
ples considered below. It seems reasonable to interpret 
ratios of successive contributions close to one as a signal 
that one needs to go beyond regular perturbative theory. 
We are getting close to this situation. For instance, the 
electro- weak corrections to — 2, give a contribution 
| p3[ of 151(4) X 10^^^ and in this calculation, the two- 
loop effects reduce the one- loop prediction by 35 percent. 
The total electroweak contribution is about one third of 
the discrepancy of 43(16) x 10^^° found by the recent 
Brookhaven experiment [Q. The problem is more se- 
rious in the case of QCD corrections. For instance, in 
the calculation n5| of the hadronic width of the Z^, the 
term of order a^is more than 60 percent of the term of 
order and contributes to one part in 1,000 to the total 
width. 

We claim that introducing large field cutoffs leads to 
significantly improved perturbative series. An important 
reference to understand the general mechanism and to 
interpret the results presented below is the well studied 
|,| integral 

/ + 00 
rf0e-(i/2)0^-A0^ . (1) 
-oo 

If we expand e"^''^ , the integrand for the order p con- 
tribution is e~(^/^)''^ (jj'^P/pl and has its maximum when 
0^ = 4p. On the other hand, the truncation of e"^"^ 
at order p is accurate provided that A^'* << p. Re- 
quiring that the peak of the integrand for the p-th order 
term is within the range of values of 4> for which the p- 
th order truncation provides an accurate approximation, 
yields the condition A << {16p)~^. One sees that the 
range of validity for A shrinks as one increases the order. 
We can avoid this problem by restricting the range of in- 
tegration in Eq. (^ to \(f)\ < 4>max- We call the truncated 
integral Z{X,(j)max)- As the order increases, the peak of 
the integrand moves across (j)max and the contribution 
is suppressed. It is easy to show that the coefficients of 
the modified series satisfy the bound \ap\ < ^/2TT(/)fP^^/pl 
and the modified series defines an entire function. How- 
ever, we are now constructing a perturbative series for a 
problem which is slightly different than the original one. 
This procedure is justified from the fact that the error is 
controlled by the inequality 



|Z(A) - Z(A,0™a,)| < 2e-^*™»^ / d</)e-(i/2)^' . (2) 

We have applied large field cutoffs to the anharmonic 
oscillator with an Hamiltonian H = p'^/2 + (jP /2 -\- \(jp . 
We use units such that the mass, the frequency and H are 
unity. The method of Ref. was used to obtain a solu- 
tion of the time-independent Schroedinger equation for 
an arbitrary value of the energy E. The eigenvalues are 
determined by using Sturm-Liouville theorem to moni- 
tor the "entrance" of the zeroes of the wave function in 
a region Q < 4> < 4>max as E increases. If ipmax is large 
enough, one obtains excellent numerical values of En for 
n not too large by requiring that the n + 1-th zero occurs 
exactly at 4>max- This numerical procedure can be con- 
verted into a perturbative expansion order by order in A. 
By taking 4>max large enough, namely 8, we were able to 
reproduce accurately the first twenty terms of the series 
for the ground state calculated by Bender and Wu |Tc| ] 
without using diagrammatic techniques. We have also 
applied large field cutoffs to Dyson's hierarchical model, 
an approximation of lattice scalar field theory where the 
renormalization group transformation can be calculated 
numerically with great accuracy We have used a 
local Lagrangian density of the Landau- Ginzburg form 
—AcfP — A0^ such that when A = 0, the mass is unity. 
The field cutoff appears in the calculation of the Fourier 
transform of the local measure necessary for the numeri- 
cal procedure. The free parameter (called c/4 in Ref. ||]) 
appearing in the kinetic term was chosen in a such a way 
that a free masslcss field scales as it would for a nearest 
neighbor model in 3 dimensions. 

We now present numerical results concerning the per- 
turbative series for Z{X, (pmax), the ground state of the 
anharmonic oscillator and the zero-momentum two point 
function of the Landau-Ginzburg model. In Fig. y, we 
compare the accuracy of regular perturbation theory with 
what is obtained with various field cutoffs. The curves for 
the modified series reach an asymptotic value on the left 
and drop on the right with the same slope as regular per- 
turbation theory but with an intercept on the coupling 
axis shifted to the right. In between the two regimes, the 
curve reaches a maximum. As the order increases, the 
maximum moves up and right making the convergence 
apparent in this region. In all cases, the modified models 
provide an accuracy that goes far beyond the envelope 
of regular perturbation theory. We have then compared 
our results with those obtained with Fade approximants. 
For the clarity of the figure, we have limited ourselves 
to order 3 and 4 but similar situations are observed for 
higher orders. At each order, we have selected the best 
approximant for the series or its Borel sum. These sums 
are obtained by dividing the l-th coefficient by r[/ -I- 1-1- 6]. 
We have followed the procedure of Ref. except that 
at the end the inverse integral transform was performed 
numerically at fixed value of A. We call this procedure 
the Pade-Borel method. 
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FIG. 2. Number of significant digits obtained with regular 
perturbation theory at order 1, 3, 5, 15 (black) and with 
(j)max = 3 (green), 2.5 (blue) and 2 (red), at order 1, 3, 11 
, as a function of A, for the three quantities described in the 
text. Even orders have high cusps and are not displayed. 
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FIG. 3. Number of significant digits obtained with field 
cutoffs given in the text, at order 3 (blue line) and 4 (red 
line) compared to the best approximants for the regular se- 
ries at order 3 (blue dots) and 4 (red dots) or the best results 
obtained with Pade-Borel method at order 3 (purple) and 4 
(green). 

The values of b were adjusted to get the best possi- 
ble result. At order 4, the best approximant is [2/2] in 
the 6 cases considered, but at order 3, the situation is 
more complicated. In summary, we have used our knowl- 
edge of the exact result to get the best possible result 
for the methods to which we compare our results. Ran- 
dom choices of approximants or of the value of h lead 
to significantly worse results. The results are shown in 
Fig. H We have used field cutoffs of 1.5 and 1 for the 
integral and 2 and 1.5 for the two other models. As the 
field cutoff decreases, the curve moves right as in Fig. |^. 
By comparing the three methods at the same order, we 
see that beyond a certain value of A, the method used 
here outperforms the two methods based on Fade ap- 
proximants within a certain range (which broadens when 
the coupling increases). 

DOUBLE-WELL 
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FIG. 4. Number of significant digits for the double-well 
at order 3 to 6 for regular perturbation (black) compared 
to series obtained with ymin = —3 and ymax = 3 (blue) or 
Vmax = 2.5 (green). As the order increases, the black curves 
reach the one-instanton contribution (red) over wider regions 
to the left while the two other sets reach the accuracy level 
obtained numerically for ymax = 3 (purple) or ymax = 2.5 
(brown) . 

In all the examples considered above, the Borel sum 
has no singularity on the positive real axis. One can in- 
troduce such singularities by adding a cubic interaction 
with an appropriate coupling. However, for all examples 
worked out, this modification can be handled properly 
with the proposed method. This is due to the fact that as 
in the previous examples, the exponentials converge uni- 
formly in a compact neighborhood of the origin, and it is 
legitimate to interchange the sum and the integral. We 
report here the case of the double-well potential in quan- 
tum mechanics as discussed in Ref. [^|. In shifted coor- 
dinates, the potential reads {l/2)y^ — gy'^ + {g^/2)y'^. By 
imposing the vanishing of the wave function at y = ±10, 
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we were able to reproduce all the significant digits of the 
first 10 coefficients for the ground state given in Table I 
of Ref. This series is not Borel-summable. We have 
then constructed a modified series by imposing the wave- 
function to vanish at ymin < and its derivative to van- 
ish at i/max > 0. If this prescription is used for numerical 
purposes, one obtains arbitrarily accurate results when 
g = {l/2ymax) and ymin negative enough. The numeri- 
cal results are shown in Fig. |^ for values of g where the 
one-instanton contribution accounts for most of the dis- 
crepancy obtained with the regular series. The modified 
series converge rapidly to the numerical value obtained 
with the corresponding ymin and ymax- It takes into ac- 
count instanton effects and significantly improves the reg- 
ular perturbative answer. Significant improvements have 
also been obtained for larger g by decreasing ymax- 

The striking resemblance among the three models ap- 
pearing in Figs. H and || suggests that, in general, the 
corrections due to the field cutoffs can be expressed as 
simple one-dimensional integrals. The perturbative ex- 
pansion of the partition function of an arbitrary lattice 
scalar field theory with a large field cutoff can be obtained 
by writing the truncated integral at each site as the in- 
tegral over the whole real axis minus the integral over 
\4>\ > 4>max- Regrouping the contributions with 0, 1, . . . 
large field contributions, we obtain the partition function 

z[j] = ce-^i:.(^)^^i:„.''fe)G(..-)^w X 

y J\'t>v\>'t>ma^ 

with Ai = (27rG(0,0))-i/2, A2 = (2G(0,0))-i, G{x,y) 
being the two-point function at A = (with no field cut- 
off) and all quantities being written in lattice spacing 
units. The dots in Eq. ^ are calculable and exponenti- 
ate in the limit of a coarse lattice where the correlations 
among the sites are small. In general, Eq. ^ can be 
interpreted in terms of Feynman diagrams. A contin- 
uum version of this expression can be obtained by using 
a dilute-gas approximation for configurations with one 
"lump" of large values. We carried the detail of this 
calculation in the case of the anharmonic oscillator using 
the classical configuration 0maa;e~'^~^°' and adapting the 
arguments of Ref. [0 . The result for the zero-th order 
correction to the ground state reads: 

/'+00 

<54°)^47r-V20L. / d<t>e-^\ (4) 

J <t>ma.ic 

This prediction fits the numerical data over a wide range 
of 4>max- We can estimate the optimal value of 4>max 
without knowing the numerical answer. The left part of 
the curves shown in Fig. || can be estimated by semi clas- 
sical methods while the right part is given by the next 
order contribution. In the case of the anharmonic oscil- 
lator, using the classical configuration mentioned above, 
we obtain the error at order n : 

\5Eo{\)\ ~ 5£;(°)e-(i/2)^^™"^ + |a„+i|A"+i (5) 



This approximate formula fits the data very well if (pmax 
is not too small and allows a good estimate of the value 
of the coupling where the accuracy peaks. 

The semi-classical calculation performed for the anhar- 
monic oscillator can be extended to other scalar theories 
with exponentially decaying two-point functions and we 
expect an exponential control of the error for these mod- 
els. In the case of lattice gauge theory, the integration 
over the fields is already reduced to a compact space. 
In the literature on lattice perturbation theory (with the 
exception of van Baal ||l^), one usually replaces / dg by 

dA], since in the continuum limit the range becomes 
infinite. We claim that this lattice artifact can be used 
to obtain a smooth truncation of the perturbative se- 
ries as in the scalar case. Different approximations need 
to be developed to solve the massless quadratic theory 
with a field cutoff. We expect that, as in the case of 
the double-well, this approach will lead to a quantitative 
understanding of the instanton effect. 
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